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Abstract 



We compute the polarization of the final-state baryon, in its rest frame, in low-energy meson-baryon 
scattering with unpolarized initial state, in Unitarized BChPT. Free parameters are determined by fitting 
total and differential cross-section data (and spin-asymmetry or polarization data if available) for pK~ , 
pK+ and pn + scattering. We also compare our results with those of leading-order BChPT. 



1 Introduction 



The study of spin phenomena in meson-baryon low-energy scattering provides stringent tests of QCD and its 
associated effective theory, Baryon Chiral Perturbation Theory (BChPT) [BUI- Because mesons are spinless, 
and at low energies can be considered nearly structureless, their scattering off baryons is the simplest process 
. from the point of view of baryon spin dynamics. As such, it is of great interest as a probe of, and may lead 

to important insights into, the structure and dynamics of baryons. 

The applicability of BChPT is limited to the near-threshold energies at which meson momenta are 
much smaller than the chiral symmetry breaking scale. At moderately higher energies, resonances and 
coupled-channels effects enter the dynamics that must be either incorporated into the theory or dynamically 
generated by it. In the three- flavor case those phenomena may be convolved with a strong-coupling regime 
• , originating in the large masses of strange hadrons. A well-known example is NK scattering, in which several 

strongly-coupled channels are open at threshold, leading to a subthreshold resonance, A(1405), and rendering 
BChPT inapplicable to those processes. Many models and techniques have been developed to overcome those 
difficulties over a period of several decades, that we cannot review here. In the specific context of BChPT, 
unitary coupled-channels techniques based on Lippmann-Schwinger or Bethe-Salpeter equations have been 
succesfully applied to the study of NK and other meson-baryon processes, even at relatively high energies 
[SI [U El [6l El [8] . A unitarization method dealing directly with the chiral effective theory T-matrix has been 
introduced in [S] in the meson sector, and extended to the baryon sector in [TOJ [TTJ [TJ] . This Unitarized 
Baryon Chiral Perturbation Theory (UBChPT) has been shown to give accurate descriptions of unpolarized 
cross-section data in NK processes [11] Q21 H31 [14] , and in Ntt scattering beyond the A resonance peak [10] . 

In this paper we consider a particular aspect of hadron spin dynamics, the production of polarized 
baryons in unpolarized meson-baryon scattering. Specifically, we compute the polarization of the final-state 
baryon in its rest frame in low-energy two-body meson-baryon scattering with unpolarized initial state, in 
UBChPT. By low energy we mean incident-meson momentum qi a b < 300 MeV. We use tree-level BChPT 
partial waves unitarized with the method of [TO] [HI [H] , and determine their free parameters by fitting total 
and differential cross-section data (and spin-asymmetry or polarization data if available) for pK~ , pK + and 
pir + scattering. We also compare our results with those of a previous leading-order BChPT calculation [T5] , 
with the aim of further probing its domain of applicability. 
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There is some unavoidable overlap with previous works (e.g., pj 111! I12j ) since we can only fit our 
calculations to the same available data as used in the previous literature. We remark, however, that although 
the unitarization method we use is the same as in [111 1121 116j , our approximation scheme is different because 
we have to include the full s- and u-channel contributions to scattering S and P partial waves, as well as 
baryon decuplet contributions, which are quantitatively important for polarization observables in the energy 
range considered here. By contrast, in order to obtain a good description of unpolarized cross sections (for, 
e.g., pK~ scattering), it is enough to use less detailed approximations. Good fits to total cross sections 
can be obtained by considering only 5-wave scattering, as done at 0(q) in pj] and at 0(q 2 ) in [14]. 
In [121 116] unpolarized differential cross sections are described in terms of S and P partial waves in a 
non-relativistic approximation in which some of the contributions mentioned above are of subleading order, 
therefore neglected. This is to be expected, since polarized observables are more sensitive to smaller partial 
waves than unpolarized ones, the latter being usually dominated by large, resonant waves. 

In the following section we present our notation and conventions, give the explicit form of the tree-level 
partial waves used throughout the paper, and very briefly discuss the unitarization method applied to those 
partial waves. In sections |3]— [5] we describe our results for pK~, pK + , and pir + scattering, resp. Detailed 
fits to scattering data are reported, and the resulting parameters applied to the computation of final-state 
polarization. In section [6] we give some final remarks. 

2 Partial waves and unitarization 

The ground-state meson and baryon octets are described by standard |17j traceless 3x3 complex matrix 
fields (f> and B, resp., with <p hermitian. We use the physical flavor basis 

ft 1 = —7= (A 1 + i\ 2 ) , /3 2 =^, /3 3 = A 3 , 

(1) 

/? 4 = -l(A^A 5 ), /3 5 =/3 4 t, /3 6 = -L(A 6 + zA 7 ), (3 7 = P 6 ^ (3 a = X 8 , 

where A a are SU(3) Gell-Mann matrices. The real matrices (3 a are not hermitian. Their hermitian conjugates 
form a basis that differs from {/3 a }a =1 only in its ordering. To distinguish field components with respect 
to each of those bases we use lower flavor indices for f3\. Thus, meson and baryon fields are decomposed 

<j> = Eb M b /V2 = £ b b /3t b /V2 and B = J2 a B a P a /V2 = E a B*fiJy/2, with <t> b = Tr (/3V) /yft, 
4> b = Tr (/3 f, 0) /V2, and similarly B a and B a . Baryon and meson states are denoted \B a (p, a)) and \M b (q)) : 
resp. , with a = ±1/2 the spin along a fixed spatial direction in the fcrmion rest frame, and p, q four-momenta. 
We always use hadron kets with an upper index, and bras with a lower one, with masses m a and fhb for 
baryons and mesons resp. Free fields couple to one-particle states as {Q\B a {x)\B c {p 1 a)) — S^u(p, a) exp(-ipx) 
and (0\(/) b (x)\M c (q)) = 5 c b exp(-iqx). 

Indices can be raised or lowered by means of the symmetric matrices e ab = Tr (f3 a f3 b ) /2 = Tr (/3 t a /3 t b ) /2 
= e ab and = 5%. In this basis the structure constants and the anticommutator constants in the fundamental 
representation are 

r,/3 b ] = 2]Tr" c /3 c , / ab c = iTr(/3t c [/3°,^]) , 

4 1 ( 2 ) 

{/3 a ,/3 fc } = 2^d Qb c /? c + | e afc , d ab c =-Tr(^ c {f3 a ,[3 b }). 

c 

Similar definitions hold for J Q b c , d a bc , etc. The constants f abc and d abc (as well as f a bc = —f abc and 

d a bc = d abc ) are totally antisymmetric and symmetric, resp. Their numerical values are different from their 
Gell-Mann-basis analogs. 

The Lagrangian of fully relativistic Baryon Chiral Perturbation Theory (BChPT) is written as a sum 
L = Cm +£mb of a purely mesonic Lagrangian Cm and a meson-baryon one Cmb- The mesonic Lagrangian 



1 0(q n ) denotes a generic quantity of chiral order n, with q a nominally small quantity such as a meson momentum or mass. 
2 We do not use summation convention for flavor indices. 
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to 0(q A ) was first obtained in [T51 US] • The meson-baryon Lagrangian Cmb has been given to 0(q 3 ) in the 
three-flavor case in [5D1 [5TJ [55] , an d hi [22] for two flavors. The tree-level amplitudes from C for meson- 
baryon scattering have been given in [TT] (see also [HHH])- We discuss here the associated partial waves, 
as used below to fit data. Some, but not all, of the S and P partial waves have been given explicitly before 
[TJJ. 

Defining the T-matrix as S — I + i(2ir) 4 5(Pf — Pi)T, the scattering amplitudes are given by T-matrix 
elements T a a r b r{s, u; a, a') = (B a ,(p', a')M b ,(q')\T\B a (p, a)M b {q)) as functions of the Mandelstam invariants 

s = (p + q) 2 , u = (p — q') 2 and the spin variables. The center-of-mass (CM) frame partial waves fg± b ,., 
corresponding to j = i ± 1/2 are defined as, 

oo 

T«\, b ,(s, u; a, ✓) = £ { ((/ + l)hf a , b , + tf t J* ' a , b ) Pe(P ■ P')X$ ■ X, 

e=o (3) 

+ * (ft+ \>« ~ h-\> b ) P e(P ■ ■ (* ■ P A P')X« } , 

with p ■ p' — cos 9c m, Pi and Pi the Legendre polynomial of order I and its derivative, and \ a , x' a ' ^~ 
component spinors for the initial and final baryon, resp. With this definition the partial waves for the 
contact interaction amplitude are, explicitly writing flavor coefficients, 

rr- ab X * b rda 

•M c ) a'b' = 2-^1 * b' d) a' ' 

d 

^ ab ab (^) 

, ab ( c ) a'b' Ar / r m a + m a <\ ab ^(c) a'b' \p'\\p\ ( r m a + m a 

Hc)o a'b' = J2 N "' N a {V* £ J , /(c)!- oV = J 2 ^V^ + 

with N a = \J p° + m a , and / the common pseudoscalar octet decay constant. For the s-channel amplitude 
we have, 

^s) ab a' b ' {d} = (Dd ba d + Ff ba d ) (Dd d bW Ff d b , a , 



a 'b' = -~F2 N a'N a (Vs ~ ma')(V^ ~ m a )Y^ 7= 

J j_i V s 



ab {d} 
(») a'b' 



m d ' (5) 



d=l 

g t- ab {d} 

t ab l \v'\\p\ , r , w r . ^ (s) a ' b ' 

where M a is the bare mass of the exchanged baryon of flavor a, and D, F are the leading-order baryon-meson 
coupling constants in SU(3) BChPT [201 [HI [22- f( c )£± an d f( s )e± vanish for all other values of £, j. For the 
w-channel partial waves we need to introduce some notations besides the corresponding flavor factor, 

:F H ab a '/ d} = (Dd\, d - Ff\^ (Dd bd a , + Ff hd a )j , 7od = m a , +m a -m d , 

2m 2 d - (u max + u min ) (6) 
lid = \™>a> + m d )(m a + m d ) , z d = , 



where the kinematic limits for u, u min < u < w max , are given as functions of s in appendix [Al The S wave is 
then writtens as, 



ab N a 'N a _ ab {d} f r- . \fs + 7od 

Hu)o a'b' - f 2 A") a'b' < V S + m d - J ld ^ ^ 

J d=1 L a' a 



(7) 



l „ V s ~ lad . V s + lod , n , v ■ 
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In this equation, and in what follows, Qi(z) denotes the Legendre function of the second kind |24j . analytic 
on the z plane cut along — 1 < z < 1 for I a nonnegative integer. In this respect, we notice that Zd > 1 
always, since octet baryons are stable under strong interactions. For the P-waves we have, 



r ab 

a'. 



1 ab {d} \\p'\\p\, r , . 7id(ys + 7od) ( n i s n , . 

d—l K 



2\p'\\p 



-jidQi(zd) 



's + 7od + 



y/s + 7od 



Aa6 
a'. 



1 I JV ' iV 

E ^^a'b^ | 2 |^/||^| (V^ ~ 10dhldQl(z d ) 



(8) 



TlrfOjXfrf) + 10d 

4\p'\\p\ N a ,N a 
Finally, for the higher u-channel partial waves, I > 1, we have 
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7irf(\/g + 7od) 
6N a ,N a 



Q2(zd) - Qo{zd) 



^s + lod. 



f ab 1 V"^ T ab Wt 1Y« ( JldQi(zd) ( AT AT ( r \ i 

I Qt+i{zd) - Qi-i{zd) ) 



x to 



/• ab 
H u )l+ a'b' 



1 r ab {d} ( - ^ / 7lriQl( z ri) Z'' at at / y x , yj + 70rf 

/ 7 lrf (^ + 7o rf )^ +i(zd) _ Q ^ ( ^ 



(9) 



2^+1 2iV -JV a 



We turn next to decuplet baryon interactions. Those have been the subject of an enormous literature 
that we cannot review here. Our starting point is the relativistic Lagrangian for A-N-n interaction from 
[25] (see also, e.g., [26[ [27J [28]). At leading chiral order the transition from two to three flavors amounts to 
inserting in the amplitudes the flavor factors for the coupling of two octets and a decuplet. The tree- level 
scattering amplitudes so obtained differ from those for Ntt scattering [10J only in their flavor coefficients. 
The CM frame partial waves for tree-level s-channel decuplet exchange are, 



r ab 
/(s,dec) a 'b' 



3 g 2 w (D + F) 2 J?* S a , b , c S a \ 



E 



s 3 / 2 m?, 
C=i c 



~\J (V^ + rn a >) 2 - to 2 , ^ (y/s + m a ) 2 - m 2 



m c 



r ab 

/(s,dec)i_ 



32 p 

\/s ((2k — l)s — 2K,\[sm a i + m\i — to 2 ,^ (j2,K — l)s — 2K\fsm a + to 2 — m 2 ^j + 

4/c(l - 2k)s 2 + (8k 2 - 4k + l)(m < + m a )s 3/2 
+ (to 2 , + to 2 — m 2 a , — m 2 a + 4k(1 — 2n)m a im a )s 

+ (-ma>ml - m 2 a ,m a + m a fhl, + m a >ml)y/s + 2(m 2 a , - m^)(m„ - m 2 ) j , 
3 g 2 (D + F) 2 ^ Sa'b'cS 



32 



E 



V (V* - m a t) 2 - ml,\J (s/s - m a ) 



2 — TO 2 X 



/s ^(2k — l)s + 2ny/sm a i + to 2 / — to 2 /^ ^(2k — l)s + 2n\fsm a + to 2 — to 2 ^ + 
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-4k(1 - 2k)s 2 + (8k 2 - An + l)(m a > + m a ) 



,3/2 



(to 2 , + m 2 — to 2 , — m 2 + 4k (1 — 2/{)m a /TO a )s 

-{-m a ,m 2 a ~ m 2 a ,m a + m a rh\, + m a 'fhl)y/s - 2(m 2 , - 777 2 ,)(to 2 - 777 2 ) 



(10) 



oh 3 g 2 p + F) 2 ^ S a , b , c S a c 

J(s,dec)i + 'ft' RzL f2 2^ 



/(s,dec) 2 - a 



64 / 2 ^ s 3/2(^_ Mc) 
Y - m a ') 2 - m 2 , y (yfs - m a ) 2 - ' 
3 g 2 w {D + F) 2 y, <Wc<S ab c 



- m h x 



64 / 2 



f lS 3 / 2 (V^+Mc) 



((Vi - 777 a ,) 2 - TO 2 ,^ ^(-y/s - 777 a ) 2 - 777^ X 



x \J{^fs + m a ,) 2 - m 2 , \j{^fs + m a ) 2 



All other partial waves f( s ,dcc)t± b a > b i with £ > 2 vanish. In (fit)]) capital Roman letters are used for flavor 

decuplet indices. Mc denotes the bare mass of the C th decuplet state. The flavor coefficients for the initial 
and final vertices in s-channel diagrams are, 



c — ^ilm{P a )jl{P b )km{Tc)ijk , Sa'b'C = ^am{Po\>)lj{Pli)m,k{Tc)ijk 



(11) 



with repeated indices i, j, . . . , summed from 1 to 3, and the matrices (3 from Tq in (fTTj) are a standard 
basis for the decuplet representation space (as given, e.g., in eq. (9) of [2"9"]). A coupling constant 510 has 
been introduced in the partial waves (fT0|) , allowing for departures from the SU(6) symmetric case <7io = 1. 

The partial waves in (|10j) were obtained by expanding the s-channel amplitude, aside from spinor factors, 
in powers of s and t and then projecting them onto the appropriate Legendre polynomials. The same 
procedure has been applied to u-channel amplitudes, related to s-channel ones by crossing, yielding partial 
waves expressed in terms of Legendre functions Qt.. The resulting expressions, however, are considerably 
lengthier than their octet baryon analogs (J?])-©, so we shall omit them for brevity. 

Tree-level BChPT is not sufficient by itself to describe three-flavor meson-baryon dynamics, in which the 
large masses of strange hadrons induce large meson-baryon couplings. This problem is especially severe in 
the NK sector, in which strong coupling effects such as subthreshold resonances render BChPT inapplicable 
even at threshold. We unitarize the tree-level amplitudes with the method first introduced in [5] in the 
meson sector and extended to meson-baryon interactions in [lOl [XT], [12j . A technically detailed explanation 
of the method can be found in those references. We shall limit ourselves here to stating the result of the 
unitarization of tree-level amplitudes. 

Given a set of coupled reaction channels \B ai M hi ) — > \B aj M bj ), we denote (fe±)ij — ft± lhl a b the 
corresponding tree-level partial wave matrix. A solution to the unitarity equation for T, resumming the 
right-hand cut in the s-plane, is given by the partial waves [Ti±)ij related to (fi±)ij by the matrix equation, 

Tt± = (I + fi± -Gy 1 -h± , (12) 

where / is an identity matrix, and G is the diagonal matrix Gij — g aibi Sij (no summation over i, j). The 
"unitarity bubbles" g ab are given by, 



f b (s) 




(13) 



with w(x, y, z) defined in appendix [XI The loop function g ab was computed in (fT3"]) in dimensional regulariza- 
tion. The subtraction constants a ab , depending on the renormalization scale fi, are taken as free parameters 
in each isospin channel. Variations in /i can be offset by a redefinition of a ab |10j . 

For numerical computations we use physical meson and baryon masses, and coupling constants D = 
0.80±0.01 and F = 0.46±0.01 [30] (see also [SI [32]). Following [TQQJ], below we set A* = 630 MeV - With 
our notation for initial and final momenta, we define the "lab" frame by p=0, the "lab"' frame by p'=0, and 
the CM frame by p + q =0. 

3 Results for pK~ scattering 

In this section we discuss our results for the reactions pK~ — > B a M b . Among the ten possible final states 
pK~ , nK , Air , E°7r°, £~7r + , £ + 7r~, An, S°?7, S ^ , E~K + , only the first six are open for initial meson 
momentum < qi a b < 300 MeV. Following [11] (see also [T6] [T2] [6]), we unitarize the tree-level partial 
waves with the method of [TO] [11] [12] taking into account all ten intermediate meson-baryon states. This is 
justified by the fact that in lowest-order BChPT those states are degenerate. Notice, however, that physical 
masses must be used in the unitarity relations, which are exact. 

In the computation of observables such as cross sections and polarizations we include the S, P1/2 and 
P3/2 partial waves arising from the contact meson-baryon interaction and s- and u-channel octet baryon 
exchange, as well as s-channel decuplet baryon exchange, in fully relativistic BChPT as discussed above. At 
the low energies considered here the contribution of D and higher partial waves is negligible. We do not 
expect u-channel decuplet baryon exchange to be significant either, as we have checked to be the case for 
several processes in this sector. 

In the isospin limit the unitarity bubble g ab contains six subtraction constants a ab (with ab — NK, £71", 
Air, £77, A77, EK), that must be fitted to data. Also extracted from data are the values of the common 
pseudoscalar octet decay constant / and the bare masses Me and Ma entering the s-channel contribution 
to the P\/2 wave, and Ms* entering the s-channel contribution to the P3/2 wave. The inclusion of baryon- 
decuplet interactions introduces two additional free parameters, the off-shell parameter Z (or equivalently, 
ft = 1/2 + Z) and the coupling 510. Whereas the latter may be made flavor-dependent, we prefer not to 
do so to avoid an excessive proliferation of parameters, and also because numerical checks suggest that 
no significant improvement of our fits is obtained by doing so. Physical masses are used in all numerical 
computations. 

In all cases we fit data up to incident momenta q\ a ^ ~300 — 320 MeV, above which D-w&ve effects start 
being important. Since fit parameters are highly correlated, we cannot vary them arbitrarily. Instead, in 
order to check the stability of, and to estimate the uncertainty in, our results we fitted the experimental 
data in several different ways, as described below. 

The data included in our fits consists of the threshold branching fractions 7, R c and R„ [33] 02], cus- 
tomarily included in pK~ coupled channel analyses [3] [Jj [11] [6] [12], total cross sections for the six open 
reactions channels, differential cross sections for pK~ — > NK, and the first two Legendre moments of the 
CM frame differential cross sections and spin asymmetries for pK~ — > £71", Air. We do not include in our fits 
data on kaonic hydrogen energy levels and their widths 35, 36 . The most recent of those experiments [36] 
has obtained precise energy and width measurements whose phenomenological consequences are discussed 

in 071 [M [M SOI US S21 El] ■ 

We present results from four different fits to data. First, we consider a non-relativistic HBChPT ap- 
proximation in which only the contact interaction contributes to the S wave, only contact interactions and 
s-channel octet baryon exchange contribute to the P1/2 wave, and only the pole term of s-channel decuplet 
baryon exchange is included in the -P3/2 wave. There is no dependence on k in this approximation, giQ is 
set to 1 and the values of the other parameters are taken from [12] . We denote this approximation as "NR" 
below. In all other fits we consider the partial waves in the fully relativistic theory as described above. 

Second, we consider a fit to the fractions 7, R c and R n , total cross-section data for the six channels open 
at giab < 300 MeV, and differential cross-section data for pK~ — * pK~ , nK . We obtain the parameters, 

/ = 106.90, a N * = -1.98, a A * = -3.55, = -1.65, a A '' = -2.30, a a ' = -5.53, a 3K = 0.20, (14) 
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with bare masses My = 1237, Ma — 1258, My* = 1320 MeV, and decuplet parameters k = 0.4 and 
gio = 1.10. We refer to this fit as "III" below. 

Third, adding to the data in fit III the Legendre moments data we obtain the parameter set 

/ = 101.44, a 1 ™ = -1.94, a AlT = -3.69, a S7r = -1.47, a Ar > = -2.11, a Sl ' = -4.04, a 3K = 5.02, (15) 

with bare masses My = 1297, Ma = 1258, Ms« = 1301 MeV, and decuplet parameters k = 0.45 and 
.9io = 1-26. We refer to this fit as "II" below. 

Lastly, we repeated the series of fits leading to II with a modified x 2 function. In order to prevent x 2 
from being dominated by the more numerous total cross section data we weighted each data set by the 
inverse number of points. With this, we obtained the slightly modified parameter set, 

/ = 102.42, a NW = -1.95, a Av = -3.88, = -1.49, a A '' = -2.13, aF v = -4.05, a 3K = 5.13, (16) 

with bare masses My, = 1127, Ma = 1258, My* = 1299 MeV, and decuplet parameters k = 0.33 and 
510 = 1.28. We refer to this fit as "I" below. 

Some remarks are in order regarding these parameter sets. First of all, we stress here that we are fitting 
data from many different experiments, which are not necessarily fully compatible among them. The values 
of / in all fits are intermediate between f v and /k, as expected |12j . The constant gio measures the coupling 
of £(1385), the only decuplet state exchanged in the s channel. In all cases we get gio ^ 1, by up to ^25%; 
these results should be taken with caution since, as per the caveat about data above, our fits are not meant 
as a determination of decuplet couplings. The range of values spanned by the subtraction constants a ab in 
fit III can be narrowed if higher x 2 values are accepted, but we prefer to reproduce the data as accurately 
as possible. Inclusion of Legendre-moment data in fits I and II leads to a noticeable dispersion in the values 
of a ab . We have found that attempting to narrow their range unacceptably worsens the resulting fits. 

Very tight fits to branching fractions and total cross sections are obtained if only those data are considered, 
but the resulting fits lead to a poor description of differential cross section data, especially Legendre moments 
for pK~ — » Stt, An. Inclusion of all differential cross section data, as in I, II, leads to slightly looser but still 
acceptable fits to total cross sections. For fit I we get the branching fractions 7 = 2.221 (2.36± 0.04), R c = 
0.650 (0.664±0.011), R n = 0.208 (0.189± 0.015), in good agreement with the experimental values [3l [34] 
quoted in parentheses. Fits II, III lead to very similar values, while NR yields a somewhat better agreement 
with data [12] - Also, as a verification we computed the T,tt mass distribution from the squared isoscalar 
amplitude for NK — > Yj-rr. We obtain a stable result for the resonance A(1405), independently of the fit we 
use the resonance peak appears at 1419 MeV with a width of 28 MeV, consistent with the results for this 
channel from the very detailed study [13) . 

Our results for total cross sections are shown in fig. [T] As expected, fit III gives slightly better results 
than I. The contributions to the S wave from s- and u-channel baryon octet exchange in those fits lead to 
a significan improvement over the results from NR. Fig. 2 shows our results for differential cross sections 
for pK~ — > pK~,nK . Although the difference between III and I is more marked in this case, with the 
latter giving a better description of data, all three curves shown in the figure provide a good description of 
data within experimental errors. The results for differential cross-section Legendre moments are significantly 
different for the three fits shown in fig. [3J As expected, parameter set I yields substantially more accurate 
results than III, NR, which were not fitted to this data. Results from fit II were omitted from these figures 
for clarity, since they are virtually identical to those from I. 

Fig. Q] shows our results for the Legendre moments Bi 2 of the CM-frame spin asymmetry. Fits I and 
II give a much better description of data than III and NR, as expected since the latter two do not include 
these data. The agreement of the former two with data is excellent for the charged modes pK~ — ► Y, ± tt t . 
For the moment B1/A0 in pK~ — > A7T° all data points have a positive central value, whereas all four fits 
yield slightly negative values. This phenomenon is seen also in fig. 8 of [7], where results very similar to our 
fit I are obtained with the same experimental data. Clearly, new data would help resolve this issue. H 

3 Also, in [7] the data for B\/Aq for pK~ — > T,°n° appear reflected through the horizontal axis. Such reflection would bring 
our fits I and If to a much better agreement with data. Still, given the large errors in these data, they have a small influence 
on fits. 
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We show our results for final baryon polarization V[ ah , in fig. [H Only results for fits I and II are shown, 
since only those give a good description of the spin-asymmetry data in fig. 21 As seen in the figure, there 
are small quantitative differences in the results for V^, from I and II for pK~ — > pK~ , E°7r°, and slightly 
larger ones in pK~ — > S + 7r _ , Att°. We take those differences as a rough estimate of the uncertainty in our 
results. 

In the process pK~ — > nK there is a global sign change in V[ ah i: at low values of gi a b V{ ah , takes small 
but negative values over the entire range of cos 9cm, which change to positive values around q\ah ~ 100 
MeV. We do not plot polarization values below ~1%, so that sign change is not visible in the figure. As 
shown there, the difference between curves I and II is larger in this case than in the processes mentioned 
above. 

A more marked sign change of V{ ah , is seen in the plots for pK ~ — » £~7r + , with small positive polarizations 
at small q^ turning negative at qi^ ~ 150-200 MeV. As is clear from the figure, this type of rapid evolution of 
7-\' ab / with energy is sensitive to the parameters and not predicted reliably within the present approximation. 
Taking into account both higher order corrections in the perturbative amplitudes and higher partial waves 
in the unitarized ones would allow us to take into account higher energy data (up to gi a b ~ 500 MeV if D 
waves are included [7j). This would result in tighter constraints on the energy dependence of amplitudes 
and, presumably, reduce the uncertainties in the extrapolation to lower energies. 

4 Results for pK + scattering 

No s-channel interactions or resonance formation are present in pK + scattering. Associated production of 
A occurs about 350 MeV above threshold in the CM frame, with the rest mass of AK at the A peak at 
<Ziab ~ 800 MeV, so it can be safely neglected in the range qi a b < 300 MeV considered here. We do take into 
account it-channel E°(1385) exchange, though its contributions can be expected to be small. We use the 
S, P\/2 and -P3/2 partial waves from contact and it-channel meson-baryon interactions, and from it-channel 
decuplet baryon exchange, at tree level as input for one-channel unitarization. 

The unitarized partial waves have four free parameters, the common pseudoscalar meson decay constant 
/, the subtraction constant a NK , and the couplings k and gio from decuplet baryon interactions, which are 
obtained from fits to differential and total cross section data [43] for gi a b < 330 MeV. Our best fit to these 
data is obtained for / = 106.04 MeV, a NK = -1.60, k = 0.4, and 310 = 0.79. This value of / is slightly 
larger than that obtained in the previous section. Almost equally good fits are obtained for values of / about 
10 MeV higher or lower, with corresponding changes in the other parameters. 

Results for differential cross sections for 145 < gi a b < 325 MeV are shown in fig. [5] We corrected our 
UBChPT partial waves for final-state Coulomb interaction with the Coulomb amplitudes of [H] (see also 
[7]). The agreement with data is excellent. Also shown in the figure is the differential cross section from 
Coulomb-uncorrected UBChPT amplitudes, and from tree-level level amplitudes including neither decuplet 
baryon nor Coulomb interactions. As expected, the contribution from it-channel baryon decuplet interactions 
is very small at these energies: setting gio=0 leads to barely visible changes in the figure. Total cross-section 
data is shown in fig. [7] 

Results for polarization are shown in fig. [3 As expected, the effect of the Coulomb interaction is to erase 
polarization, and it is stronger at lower qi^ and in the forward direction. Also shown in the figure are the 
UBChPT results without baryon-decuplet contributions, and the result from leading-order BChPT (without 
both decuplet baryon and Coulomb interactions) for VI ^ [15] . As seen in the figure the contribution to 
V[ ab , from £°(1385) exchange is small, about ~ 10% of the polarization value at the higher end of the 
energy range in the figure and smaller at lower energies. The agreement between Coulomb-uncorrected 
UBChPT without decuplet baryons and leading-order BChPT results is remarkably good up to q\ a b ~ 300 
MeV, where higher-order correction effects become apparent. The fact that the UBChPT polarization is well 
approximated by the leading-order one at low energies shows that the former scales with / approximately 
as l// 2 , as the latter does exactly. 
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5 Results for Nn scattering 



Ntt scattering at low energies, (ft a b < 100 MeV, is well described by one-loop HBChPT [JS]. Inclusion 
of A degrees of freedom extends the applicability of HBChPT up to <7i a b < 300 MeV, as shown in 46J. 
Unitarization, and inclusion of the A(1440) and p meson, allows further extension to gi a b < 400 MeV 10J. 
In the context of fully relativistic BChPT the 7V7rA interaction has been studied at one-loop level in [47] , 
No detailed phenomenological analyses of Ntt processes such as [45] [46] have been given in one-loop BChPT 
yet. In this section we apply the unitarized tree- level BChPT partial waves given above to Ntt scattering. 
At near-threshold energies, tree-level UBChPT is necessarily less precise than a full one-loop treatment. As 
discussed below, however, it yields a very good description of Ntt data up to <7i a b < 350 MeV, with many 
fewer free parameters. 

For simplicity we restrict ourselves here to pn + scattering. Besides the elastic process, also pn + — > 
Tr r K + should in principle be considered. Unitarity corrections arising from the couplings of those channels, 
separated by more than 600 MeV in the CM frame, are expected to be negligible at the energies considered 
here. Indeed, if that were not so the applicability of BChPT to Ntt scattering would be called into question. 
We therefore apply a one-channel unitarization procedure to the elastic pir + amplitude, analogous to the 
treatment of pK + scattering above. 

We unitarize the tree-level partial waves 5, -P1/2, -P3/2 computed from contact, s- and u-channel meson- 
baryon diagrams, and s- and it-channel A exchange. The free parameters are determined from a fit to 
differential cross-section and analyzing power data in the region 150 < < 330 MeV, around and below 
the A resonance. A series of fits was performed in order to check the stability of our results, with only small 
differences found among them. The parameter k, in particular, was found to be around 0.4 in all cases, with 
small fluctuations [10]. Using only differential cross section data, the best fit is found for 

/ = 78.81 MeV , a N * = -1.54 , k = 0.38 , g w = 0.82 , (17) 

with a bare mass Ma = 1285 MeV. We refer to this fit as "I" below. For elastic processes lab-frame 
analyzing power is equal to lab'-frame polarization, as can be easily checked (see, however, [IB] for an 
exhaustive discussion of polarization observables) . Including analyzing power data in our fits we find, 

/ = 78.39 MeV , a Nn = -1.77 , k = 0.37 , g w = 0.81 , (18) 

with a bare mass Ma = 1299.5 MeV. We refer to this fit as "II". The values for / obtained in these fits 
is lower than the one from fits to pK ± data above, which may be explained by the fact that only Ntt 
interactions are present in this case. With respect to the value of gio the same caveats as in sect. [3] apply 
here. The A resonance peak in the P3/2 wave is located at ^1226 MeV, with a Breit-Wigner width of ^110 
MeV, in both cases I and II. 

The resulting differential cross sections are plotted in fig. [9] for gi a b < 350 MeV. Only a representative 
sample of data is shown in the figure, including only the energies with the largest number of data points, 
since the entire set comprises twenty different values of gi a b- The agreement with data is excellent, given 
the small experimental errors especially in the most recent experiments [33] ■ For these data the difference 
between I and II is barely visible in the figure. Polarization results are shown in fig. [TO] to be in very good 
agreement with dataQ The difference between the two fits is noticeable here at peak polarizations. Aside 
from that small difference, the result I is remarkable since polarization data were not included in that fit. 

If we neglect decuplet baryons, the results from UBChPT agree quite precisely with those of leading 
order BChPT [15] up to gi a b ~ 300 MeV, similarly to the case of pK + scattering discussed above. Inclusion 
of A degrees of freedom leads to large effects. A comparison of results with and without A was given in 
[15] . where only the pole term in the -P3/2 wave was taken into account. The polarizations shown in fig. [11] 
are substantially larger than those in |15j . The latter are reproduced in the present approach if, firstly, we 
increase the value of / (here extracted from data) to the larger nominal value used in [15j . And, second, 
if we restrict the A-exchange contribution to s-channel P3/2 wave as done in that ref. The contribution of 
s-channel A exchange to the S partial wave, in particular, is responsible for the enhancement of polarization 

4 Notice that we define polarization in lab frame with respect to q Ap' oc —qAq' |15) . so a sign-flip has been applied to the 
experimental data in fig. 1101 
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with respect to |15j not accounted for by the change in /. As shown in fig. 1111 the contribution of u-channel 
A exchange to V{ ab , is comparatively small, except near the backward direction at the highest energies 
plotted there. Together with the results of [15], those presented in this figure give a complete picture of 
A-resonance contributions to nucleon polarization at low energies. 

6 Final remarks 

In this paper we computed the polarization of the final-state baryon in its rest frame in low-energy meson- 
baryon scattering with unpolarized initial state, in the unitarization framework of [101 111! 112] applied to 
fully relativistic leading-order BChPT. 

For pK~ scattering very good agreement is found between UBChPT and unpolarized cross section data, 
confirming (and to some extent, improving) several previous analyses [11 ( 1121 , [16]. The agreement with CM- 
frame spin-asymmetry data is also good, though somewhat less certain, as seen in fig. [4] We attribute this 
to the scarcity of such data and the rather large errors in some of the existing measurements. On the other 
hand, it is precisely spin observables that put the tightest constraints on fit parameters, since the clearest 
discrimination among our different fits is provided by fig. 2] 

The results in fig. [5] show large values of polarization for pK~ — > Y,tt, Att even at relatively low energies 
(l\ab ~ 150 MeV). It is interesting also that sizable polarizations are obtained for the elastic and quasielastic 
processes pK~ — > pK~ , nK° at the upper end of the energy range considered in fig. [5J The stability of 
the results, as indicated by the two curves in fig. [5J varies widely for the different processes. Whereas the 
results for pK~, E°7r°, £ + 7r~, Att final states seem reliable, those for nK are less stable, and for £~7r + 
largely inconclusive. New, more precise data on unpolarized differential cross sections at low energies would 
certainly be most helpful to reduce uncertainties. Alternatively, as mentioned in sect. [3] extending the 
theoretical computation to fit the more abundant higher-energy data (along the lines of [7]) could lead to a 
more tightly constrained extrapolation to lower energies. 

For pK + scattering excellent agreement between our UBChPT results and differential cross-section data 
is shown in fig. [SI Significant (> 10%) values of polarization are found at the upper end of the energy 
range of fig. [8] even after correcting for Coulomb interaction. The agreement of the UBChPT result for 
polarization with that of leading-order BChPT [15] in fig. [5] is quite remarkable for q\ ab < 300 MeV. Since 
those calculations were carried out by very different methods, their agreement provides a cross-check for 
both. Experimental data at these energies would obviously be most desirable. 

In the case of ptt + scattering also excellent agreement is found with differential cross-section data in the 
A peak region, as (partially) shown in fig. [9] Such agreement is remarkable given the high precision of some 
of the data sets used, and the fact that our results are based on unitarized tree-level amplitudes. More 
importantly, our results for polarization are in very good agreement with experimental data, as seen in fig. 
Ill[ which provides a non-trivial validation to our calculations. 

Interestingly, a comparison of fig. [TT] with the results of |15j shows that the non-resonant S- wave contri- 
bution from s-channel A exchange has a strong influence on polarization, suggesting that precise polarization 
data could discriminate among different A-A-7T interaction models. As also shown in the figure, the contri- 
bution to polarization from it-channel A exchange, though significant, is relatively small. 

We hope that the results presented here may motivate new measurements at low energies, or re-analyses 
of existing data, leading to new and more precise experimental information on polarization observables 
especially in the strange sector. 
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A Kinematics 



In this appendix we gather some kinematical definitions used throughout the paper. We introduce the 
notation 

u{x,y,z)= {x 2 + y 2 + z 2 - 2xy - 2xz - 2yz)? = (x - {y/y + VI) 2 ) 3 (x - (y/y - ^fz) 2 )^ . (19) 

The function u appears frequently in relativistic kinematics (e.g., in the center of mass frame \p\ = 
Lu(s,m 2 a ,m 2 ) I \2yfs~)). The Mandelstam invariants for the process \B a (p, a)M b (q)) — > \B a ' (p', a')M h ' (q')) 
are 

■s = (p + q) 2 = (p' + q') 2 , t=(p-p') 2 = (q-q') 2 , U = (p - q'f = (p' - qf , (20) 
with s + t + u = to 2 + to 2 , + m 2 + m 2 , . The physical region for the process is defined by the inequalities 



s t h < s , i min < t < i max , w min < u < u max , (21) 



where, 



s th = max 



{(m a + TO h ) 2 , (m a > + m b ,) 2 ) 



Umax 



( 2 / 2 i 2 i ~2 i ~2\ i / 2 ~2\/ 2 ~2\\ 

~2s V ~ s{jTla m °' + m b + m b') + (m Q - m b )(m a , - m v ) J 



2s 

±—uj(s,m 2 a ,ml)w(s,m 2 a ,,m 2 b ,), ( 2 2) 

~ (s 2 - S(TO 2 + TO 2 , + TO 2 + TO 2 ) - (TO 2 - TO 2 )(TO 2 , - TO 2 ,)) 



2s 

|1 c 2 ~ 2\ / 2 ~ 2 \ 
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Figure 1: Total cross sections for pK scattering. Solid lines: results from fit I, dot-dashed lines: fit III, 
dotted lines: fit NR. See text for explanation of fits. Data from [50 | I5T 1 152 } f53 j 154 } [55] . 
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Figure 2: Differential cross sections for pK scattering. Solid, dot-dashed and dotted lines as in fig. [TJ Data 
from 1541. 
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Figure 3: Legendre moments Ai t 2 for pK scattering. Solid, dot-dashed and dotted lines as in fig.Q] Data 
from [501 [55]. 
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Figure 4: Legendre moments -61,2 for pK scattering. Solid, dot-dashed and dotted lines as in fig. [T] Dashed 
lines: fit II. Data from [S01 [55]. 
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Figure 5: UBChPT results for the polarization of the final baryon in its rest frame in pK~ unpolarized 
scattering, as a function of center-of-mass scattering angle. Solid lines: parameter set I, dashed lines: 
parameter set II (see text for an explanation of parameter sets). Curves (1) — (6) (not all shown in some 
plots for clarity): q uh = 50, 100, 150, 200, 250, 300 MeV, resp. 
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Figure 6: Differential cross section for pK + scattering. Solid lines: UBChPT results corrected for Coulomb 
interaction, dotted lines: uncorrected UBChPT results, dot-dashed lines: tree-level BChPT results. Data 
from [43] . 
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Figure 7: Total cross section for pK + scattering. Solid line: UBChPT result. Data corrected for Coulomb 
interaction from [43j . 
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Figure 8: UBChPT results for the polarization of the final baryon in its rest frame in pK + unpolarized 
scattering, as a function of center-of-mass frame scattering angle. Left panel: complete UBChPT including 
baryon decuplet interactions. Right panel: UBChPT with baryon decuplet interactions switched off. Solid 
lines: UBChPT with Coulomb interaction, dot-dashed lines: UBChPT without Coulomb corrections, dotted 
lines: leading-order BChPT calculation without both baryon decuplet and Coulomb interactions [TS]. Curves 
(1)— (5): q lah = 150, 200, 250, 300, 350 MeV, resp. 
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Figure 10: UBChPT results for the polarization of the final baryon in its rest frame in pir + unpolarized 
scattering, as a function of center-of-mass frame scattering angle. Dashed lines: fit I, solid lines: fit II. Data 
from [59l [60l [61] , after sign flipping. (See main text for explanation of fits and sign convention.) For data 
from [61], statistical and systematic errors added in quadrature. 
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Figure 11: UBChPT results for the polarization of the final baryon in its rest frame in p-K + unpolarized 
scattering, as a function of center-of-mass scattering angle. Solid lines: full computation with both s- and 
u-channcl A-exchange diagrams, with parameters from fit I. Dashed lines: same computation, without 
ii-channel A-exchange diagram. Curves (1) — (6): q\ a b = 125, 150, 200, 250, 300, 350 MeV, resp. 
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